Spin gap behaviour in a 2-leg spin-ladder BiCu2P06 
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Abstract 

We present magnetic suscceptibility and heat capacity data on a new 5 = 1/2 two-leg spin ladder 
compound BiC^POg. From our susceptibility analysis, we find that the leg coupling Ji/fcg is ~ 80 
K and the ratio of the rung to leg coupling J2/J1 ~ 0.9. We present the magnetic contribution to 
the heat capacity of a two-leg ladder. The spin-gap A/ks = 34 K obtained from the heat capacity 
agrees very well with that obtained from the magnetic susceptibility. Significant inter-ladder 
coupling is suggested from the susceptibility analysis. The hopping integrals determined using 
iVth order muffin tin orbital (NMTO) based downfolding method lead to ratios of various exchange 
couplings in agreement with our experimental data. Based on our band structure analysis, we find 
the inter-ladder coupling in the froplane J3 to be about 0.75Ji placing the compound presumably 
close to the quantum critical limit. 

PACS numbers: 75.10.Pq, 75.40.Cx,71.20.-b 
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Introduction: Following the discovery of high-temperature superconductivity (HTSC) in 
the cuprates [l|, there has been an increased focus on the properties of low-D antiferro- 
magnetic systems. This is due to the innate exotic properties of these magnetic systems 
themselves and their supposed connection with HTSC. Significant work has taken place 
recently elucidating the properties of S — 1/2 and S = 1 Heisenberg chains and their re- 
sponse to impurity substitutions. Whereas quantum fluctuations prevent long-range order 
(LRO) in ID Heisenberg systems, 3D systems exhibit conventional LRO. On the other 
hand, in 2D systems where the strength of magnetic interactions and quantum fluctuations 
can be comparable, one might expect competing ground states and a quantum critical point 
separating them. Spin-ladders serve as a bridge between one-dimensional (ID) and two- 
dimensional (2D) magnetic systems and it is believed that an improved understanding of 
spin-ladders will lead to a better understanding of magnetism in the 2D systems. A major 
step was taken in this direction with the prediction of spin-gaps in even-leg ladders and 
their absence in odd-leg ladders. [2| followed by experimental verification in SrC^Os (2-leg 
ladder) and Sr 2 Cu30s (3-leg ladder).^ However, in spite of the large experimental effort, 
only a small number of gapped ladders have been synthesized and studied. Of these, only 
two (LaCu0 2 5 and Sr 14 Cu 24 4 i)could be doped significantly with holes of which only the 
latter becomes superconducting |4|. Some other compounds which have been investigated 
are (C 5 H 12 N) 2 CuBr 4 (Ref. fl), Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 (Ref. fl), Cu 2 (C 5 H 12 N 2 ) 2 Br 4 (Ref. Q) 
which have substantially smaller spin-gaps. There is continued effort to synthesize and study 
new low-D systems since they provide a rare opportunity to elucidate the significance of low- 
dimensionality, spin-gap, etc. to HTSC as also allow to examine impurity /doping effects in 
a strongly correlated cuprate. 

In this paper, we report on the preparation and properties of a new cuprate which, we 
demonstrate, can be modeled as a two-leg ladder system with significant inter-ladder cou- 
pling in the fee-plane and negligible interplanar coupling. The spin-gap, as determined from 
our susceptibility and heat capacity measurements is about 34 K while the intraladder leg 
coupling is about 80 K. Our electronic structure calculations within the framework of iVth 
order muffin tin orbital (NMTO) downfolding method 8J yield hopping integrals between 
various Cu atoms. Using the NMTO downfolding method, we calculate the Wannier-like 
effective orbitals which illustrate the shape and extent of the active Cu orbitals and therefore 
indicate the exchange pathways which lead to the ladder topology. From a practical stand- 
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point, the estimated J/k B ~ 80 K provides a unique opportunity to examine the excitations 
of the coupled ladder system at temperatures ranging from well above J/k B to well below 
J/k B . Impurity substitutions will then allow to probe the nature of magnetic effects thus 
induced, in a wide temperature range. 

Crystal structure and measurements: Our measurements are on single phase, polycrys- 
talline BiC^POe samples (space group Prima with lattice parameters a = 11.776 A, 
b = 5.1776 A, c = 7.7903 A). 

A schematic diagram of the structure is shown in Fig. 1. The unit cell contains four 
formula units, with two inequivalent Cu (Cul and Cu2) sites and four inequivalent O (01- 
04) sites. The characteristic feature of the structure are CuOs distorted square pyramids, 
with Cu 2+ ion at the centre of the five fold oxygen co-ordination. Two such pyramids 
share an edge formed from a pair of basal oxygens (02) to give rise to a Cu dimer with an 
intradimer distance of 2.8 A. Along the 6-axis, each dimer connects two others by its four 
01 corners resulting in a zigzag double chain (ladder) running along the 6-axis (see Fig. 
1). The interdimer cohesion is further strenghtened by PO4 tetrahedra that connect two 
consecutive dimers by 02 corners. The Bi ions are positioned between two ladders. The 
Cu-O-Cu angle along the leg is about 112° and that along the rung is about 92°. In Fig. 
1, various exchange couplings (J%, J2, etc.) and hopping integrals (t±, £2, etc.) have been 
indicated. 

Our results of the susceptibility Xmeas (magnetisation M divided by applied field H) as 
a function of temperature T using a vibrating sample magnetometer (VSM) of a Physical 
Property Measurement System (PPMS) from Quantum Design are plotted in Fig. 2. As 
seen, Xmeas has a broad maximum around 57 K (indicative of a low-D magnetic system) 
below which it drops rapidly (suggestive of a spin-gap). A very weak low-temperature 
upturn is seen below 6.5 K, likely due to extrinsic paramagnetic impurities and/or natural 
chain breaks in our polycrystalline sample. We now analyse these data quantitatively. 

An analytical solution for the spin-susceptibility of two-leg ladders in the full T-range 
is not known. However, Johnston, based on extensive quantum Monte Carlo (QMC) 

Q 

simulations |9J] has proposed an equation which accurately reproduces the QMC-determined 
susceptibilities at discrete temperatures. This equation (not reproduced here since it is 
unwieldy) is useful for determining the exchange couplings by fitting the measured suscep- 
tibility data and has been used to analyse such data in the two-leg ladder SrCu203. We 



3 



then fit (dashed line in Fig. 2) Xmeas to Xo + C/(T — 9) + rnxi a dder{T) where the fitting 
parameters are Xo, C, 9, J2/J1, Ji, and m. Here XiadderiT) is the x °f isolated ladders 
as given by Johnston. 9j In the absence of a generic fitting function which can take into 
account arbitrary inter-ladder interactions, we attempt to do so using the parameter m. 
With m as a variable, the obtained parameters are Xo — (4.4 ± 0.1) x 10~ 4 cm 3 /mol Cu, C 
= (3.0 ± 0.2) x 10- 4 cm 3 K/mol Cu, 9 ~ K, J 2 /J x = 0.87 ± 0.05, J x jk B = (80 ± 2) K, and 
m = 0.41±0.02. The value of the spin-gap usingfl A/J x = 0.4030( J)+0.0989(^) 3 is about 
34 K. The Curie constant corresponds to less than 0.1 % of isolated 5=1/2 impurities. 
This value is comparable to typical parasitic Curie terms found in single crystals, indicating 
the very high quality of our samples. Since the core-diamagnetic susceptibility [101] Xcore is 
—0.6 x 10~ 4 cm 3 /mol, Xo — Xcore yields the Van Vleck susceptibility xvv = 5 x 10 _4 cm 3 /mole 
which is somewhat higher than xvv of other cuprates. We show in Fig. 2 the curve for 
isolated 2-leg ladders (with Ji/ks = 80 K). We also show the simulated curve for a uni- 
form, 2D S — 1/2 HAF with J/ks = 80 K where the high-T behavior is generated using 
the series expansion given by Rushbrooke and Wood.|ll| Also, Johnston 12j] parametrized 
the low-T (^f- < l) simulations of Takahashi 13] and Makivic and Ding, b4| which we use. 
The experimental data are lower than both, the 2D HAF curve and the isolated ladder 
susceptibility. This new behavior points to the importance of a next-near-neighbor (NNN) 
interaction along the leg (which might be expected due to the zig-zag nature of the leg) which 
might be frustrating and might even enhance the spin-gap. In a latter section, based on our 
band-structure calculations, we actually find significant NNN as also inter-ladder couplings. 
The absence of LRO inspite of these deviations from the isolated ladder picture should mo- 
tivate the theorists to refine their models of such systems. In the inset of Fig. 2, we plot the 
normalised susceptibility x*(T) = x sp in(T)Ji/(Ng 2 /j, 2 B ) (where Xspin(T) Xmeas Xo 

C/T) 

as a function of fc^T/Ji. We find ^*> max ( i.e., x* at the broad maximum) to be about 0.05 
which is lower than the expected value for isolated ladders of about 0.12. 

To further confirm the spin-gap nature of BiCu 2 P0 6 , we did heat capacity C p mea- 
surements. Since the lattice C p dominates the data, it has so far not been possible to 
experimentally determine the magnetic contribution to C p in any spin-ladder compound 
unambiguously. In the present case, we are fortunate to have a non-magnetic analog of 
BiCu2P06 in BiZn 2 P06- We have then determined the magnetic heat capacity Cm of 
BiCu 2 P0 6 by subtracting the measured C p of BiZn 2 P0 6 from that of BiCu 2 P0 6 (see Fig. 3 
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inset). Ours are the first reported Cm data in a spin-ladder compound. The data are fit to 



^C M (T) = lNk B (A) 1/2 (_A_) 3/2 



1 + ^ + 0.75 



A 1 \k B T J 

k B 



exp(^A) shown by the solid 



line (Fig. 3 inset). From the fit, the spin gap — ~ 34 K, in excellent agreement with our 



susceptibility results. 

First principles study: The local density approximation-density functional theory (LDA- 
DFT) band structure for BiCu2P06 is calculated using the linearized-muffin-tin-orbital 
(LMTO) method based on the Stuttgart TB-LMTO-47 code. [l5] The key feature of the 
non spin-polarized electronic structure presented in Fig. 4 are eight bands crossing the 
Fermi level which are well separated from rest of the bands. These bands are predominantly 
derived from the antibonding linear combination of Cu d x 2_ y 2 and basal O p a states in the 
local reference frame where the z-axis is along the shortest Cu-0 bond while the x and y 
axes point along the basal oxygens 01 and 02. The band structure is 2D with practically 
no dispersion perpendicular to the ladder plane (along FX). The eight band complex is 
half-filled and metallic as expected in LDA. It lies above the other occupied Cu-d, O-p, and 
Bi-s character dominated bands. The P (s,p) and Bi (p) derived states are unoccupied and 
lie above the Fermi level, with the Bi-p states having non-negligible admixture with the 
conduction bands. This admixture of the conduction band with Bi-p states is important 
in mediating the Cu-Cu inter-ladder exchange coupling. Starting from such a density func- 
tional input we construct a low-energy model Hamiltonian using the NMTO downfolding 
technique. This method [8j extracts energy selective Wannier-like effective orbitals by inte- 
grating out high energy degrees of freedom. The few orbital Hamiltonian is then constructed 
in the basis of these Wannier-like effective orbitals. Here, we shall retain only Cu d x 2_ y 2 
orbital in the basis and downfold the rest. The effective Cu d x 2_ y 2 muffin tin orbitals gen- 
erated in the process will be renormalised to contain in their tail other Cu-cf, O-p, Bi, and 
P states with weights proportional to the admixture of these states with Cu d x 2_ y 2. Fourier 
transform in the downfolded Cu d x 2_ y 2 basis gives the desired tight-binding Hamiltonian 
7i = tij ( c j c « + c l c jj i n terms of the dominant Cu-Cu hopping integrals t; L j. This tight 

binding Hamiltonian will serve as the single electron part of the many-body Hubbard model 



relevant for this system and can be mapped to an extended Heisenberg model with the 

if 2 

■i 

J eff 

onsite Coulomb interaction. The various hoppings are displayed in Table 1 and indicated 



exchange couplings related to the LDA hoppings by = jf 3 - where U e ff is the screened 



in Fig. 1. The intra-dimer (rung) hopping proceeds mainly via the edge sharing oxygens 



5 



while the inter-dimer interaction (leg hopping) proceeds via the corner sharing oxygens with 
support from the PO4 complex. As anticipated in the experiments, we do indeed find that 
the ratio of the rung hopping to the leg hopping J 2 / J\ ~ 1. We find that the NNN coupling 
along the leg J4 is about 0.3 J\. Depending on the relative sign of this interaction with 
respect to that of J\ one might get significant frustration effects which should also have a 
bearing on the ground state of the system. We also find an appreciable coupling between the 
ladders (J3/J1 ~ 0.75) mediated primarily by the unoccupied Bi— p states. Our conclusion 
is further supported by the plot of the corresponding Cu d x 2_ y 2 Wannier function in Fig. 
5. We find that each Cu d x 2_ y 2 orbital in the unit cell forms strong pdcr antibonding with 
the neighboring 0-p x /0-p y orbitals resulting in the conduction band complex. The Cu ions 
strongly couple along the leg as well as the rung of the ladder confirming the hoppings in 
either direction should be comparable. The 0-p x /0-p y tails bend towards the Bi atom, indi- 
cating the importance of the hybridization effect from the Bi cations and therefore enhances 
the Cu-Cu interladder exchange interaction (see table 1). 

Conclusion: In conclusion, we have presented a new S — 1/2 2-leg ladder BiC^POg. 
From our x and Cm data we obtain a spin-gap Ajks ~ 34 K and a leg coupling Ji/ks ~ 80 
K. From our first principles LDA-DFT calculations, we find J2/J1 ~ 1 and a significant 
inter-ladder interaction in the corrugated fee-plane (J3/J1 ~ 0.74). Considering that the 
uniform S — 1/2 2D AF system has an ordered ground state, we feel that the strong inter- 
ladder interaction in BiCu 2 P0 6 places it close to a quantum critical point. The moderate 
value of the A/ks in BiCu2P06 will allow one to explore the magnetic properties in a large 
T range, well below and well above the gap temperature, enabling a comparison with and 
refinement of theoretical models. We feel that there might still be unanticipated features 
in the physics of low-D magnets and we expect our work to motivate others to carry out 
numerical simulations and explore the phase-diagram of coupled 2-leg ladders in the presence 
of NNN couplings along the leg. We are presently considering doping/substitutions in this 
2-leg ladder which might be able to tune the inter-ladder exchange and effect a quantum 
phase transition. 

We thank the Indo-French Center for the Promotion of Advanced Research for financial 
support. 
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Table Caption 

Table 1 Hopping parameters (t n ) between various Cu's are tabulated along with the 
corresponding Cu-Cu distances. The hopping paths are indicated in Fig. 1. 
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Hopping path 


Cu-Cu distance 
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in meV 
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\ v n / v i } 


lee (ti) 


3.22 
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runf? (to} 
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inter-ladder (t^) 
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0.74 


NNN in leg (t 4 ) 


5.18 


91 


0.34 


diagonal (£5,4) 


4.43 


30 


0.04 


diagonal (is a) 


5.81 


26 


0.03 



Figure captions 



Fig. 1 (Color Online) A schematic of the BiCu 2 POg crystal structure is shown. It 
can be seen that 2-leg ladders run along the crystallographic 6-direction. The 2-leg ladder 
is separately shown for clarity. Also shown are the various significant hopping parame- 
ters/exchange couplings between Cu atoms. 

Fig. 2 (Color Online) Magnetic susceptibility (Xmeas = M/H) vs. temperature T for 
BiCu 2 POg in an applied field of 5 kG. The open circles represent the raw data and the 
dashed line is a fit (see text). Also shown are simulated curves for the isolated ladder (dark 
gray line) and for the 2D HAF (gray line). The inset shows the dependence of x* 011 ksT/ J\ 
(see text). 

Fig. 3 (Color Online) The measured heat capacity as a function of T for BiCu 2 P0 6 and 
BiZn 2 P0 6 . Inset: the magnetic specific heat of BiCu 2 P0 6 along with a fit (see text). 

Fig. 4 (Color Online) LDA band dispersion of BiCu 2 P0 6 along various symmetry direc- 
tions. 

Fig. 5 (Color Online) Effective Cul d x 2_ y 2 orbital with lobes of opposite signs colored as 
black and white. The d x 2_ y 2 orbital is defined with the choice of local coordinate system as 
discussed in the text (height of the isosurface = ±0.09). The spheres represent the ions. 
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